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Abstract 

We discuss the theory of spin waves in non-degenerate ultra-cold gases, and compare 
various methods which can be used to obtain appropriate kinetic equations. We then study 
non-hydrodynamic situations, where the amplitude of spin waves is sufficiently large to 
bring the system far from local equilibrium. The full position and momentum dependence 
of the distribution function must then be retained. 

In the first part of the article, we compare two general methods which can be used 
to derive a kinetic equation for a dilute gas of atoms (bosons or fcrmions) with two 
internal states (treated as a pseudo-spin 1/2). The collisional methods are in the spirit 
of Boltzmann's original derivation of his kinetic equation where, at each point of space, 
the effects of all sorts of possible binary collisions are added. We discuss two different 
versions of collisional methods, the Yvon-Snider approach and the S matrix approach. 
The second method uses the notion of mean field, which modifies the drift term of the 
kinetic equation, in the line of the Landau theory of transport in quantum liquids. For 
a dilute cold gas, it turns out that all these derivations lead to the same drift terms in 
the transport equation, but differ in the precise expression of the collision integral and in 
higher order gradient terms. 

In the second part of the article, the kinetic equation is applied to spin waves (or 
internal conversion) in trapped ultra-cold gases. Numerical simulations are used to illus- 
trate the strongly non-hydrodynamic character of the spin waves recently observed with 
trapped ^^Rb atoms. The decay of the phenomenon, which takes place when the system 
relaxes back towards equilibrium, is also discussed, with a short comment on decoherence. 

In two appendices we calculate the Wigner transform of the interaction term in the S 
matrix method, to first order in gradients; appendix I treats the case of spin-independent 
interactions, appendix II that of spin-dependent interactions. 

1 Introduction 

Recent experiments ^ |2] have renewed the interest in spin waves in dilute quantum gases, 
in conditions where they had not been observed before. A first novehy is that these waves 
were neither purely nuclear nor electronic, but involved two hyperfine atomic levels where the 
nuclear and electronic spin are coupled. This difi'erence is actually minor: the exchange effets 
from which the waves originate are unaffected by the precise nature of the levels - in other 
words one can, without any loss of generality, assimilate any pair of atomic states to up and 
down states of a fictitious spin. The second, more important, difference is that the waves were 
observed with very large amplitudes, leading to an almost complete apparent segregation of 
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the atoms in each internal state, and involving situations where the gas is very far from local 
equilibrium; we study such situations in the present article. 

In degenerate liquids, the existence of spin waves has been known for many years j^lHj. In 
dilute gases, which are non-degenerate, it took about twenty more years jSHH] to realize that 
they should also sustain similar waves, for both bosonic and fermionic systems. This came 
as a surprise to some since, before, these waves had been mostly associated with the Landau 
formalism for Fermi liquids, and with the notion of "molecular field" resulting from exchange 
and interactions with nearest neighbors. In a liquid, a test particle remains trapped in a sort 
of local cavity and interacts constantly with several neighbors, so that the notion of mean field 
emerges naturally from the averaging over the effect of several neighbors. By contrast, in a gas, 
the motion of a particle is often described as free flights along straight lines, interrupted by 
short collisional processes, so that the physics involved is clearly very different. It nevertheless 
turns out that, when averaged over all possible collisions at each point of space, the final result 
of the exchange interactions in a gas and a liquid are very similar (see, for example, Ref. [Jj 
for a general discussion). Soon after these theoretical predictions, experimental observations 
demonstrated the existence of spin waves in spin polarized Hydrogen gas 8 , '^He gas ^Q,, and 
dilute ^He-'^He solutions pi7] . 

The physical difference between gases and liquids has its counterpart in the different the- 
oretical approaches used to derive kinetic equations. Historically, the first derivation was that 
of Boltzmann, with his Stosszahlansatz (or molecular chaos Ansatz): collisions are considered 
as point processes (with no duration and no spatial extent) taking place between particles 
which are completely uncorrelated before collision. As a consequence, one can add at each 
point of space the effect of all possible binary collisions between uncorrelated particles, as if 
one was adding the effect of many "beam to beam collision experiments" in atomic physics. 
This provides the famous Boltzmann equation, with a relatively complicated collision term 
on the right hand side, while the drift term on the left hand side corresponds merely to a 
completely free flight of particles between collisions. In the approach introduced later by 
Landau for the study of degenerate Fermi liquids JJ, quasi-particles are never free, since 
their motion is constantly guided by a mean field created by the neighbor particles. Therefore 
interactions modify the drift term while, on the other hand, the collision term is generally 
treated phenomenologically, by a simple relaxation approximation. Both point of view have 
their advantages. Boltzmann's point of view does not require the use of a pseudo-potential 
and therefore allows a more microscopic treatment of collisions, for instance a more precise 
inclusion of lateral collisions including their full angular end energy dependence. On the other 
hand, it does not encompass quantum degenerate systems, does not introduce the powerful 
notion of quasi-particles, and generally speaking remains limited to dilute interacting gases. 

The connection between the two points of view is provided by collisions in the forward 
direction. In two body collision theory, it is well known that wave interference in the forward 
direction is responsible for particle absorption (optical theorem). In addition, more physics 
may be involved in this forward interference effect: for spinless particles, retardation effects 
are introduced (see for instance ^,); for particles with spin, in addition, identical spin rota- 
tion effects (ISRE) due to particle indistinguishability Ej also take place. In many-body 
transport theory, both these effects have their equivalents. The cumulated retardation effects 
during collisions give rise to an average force, or equivalently to a scalar mean field (see for 
instance ^Ij and references therein); in the kinetic equation, this corresponds to terms which 
are quadratic in the distribution function, as the collision term, but contain gradients (spatial 
gradients of the density as well as momentum gradients of the particle distribution) ; in other 



2 



words, one gets zero sound Landau type mean field terms which are naturally grouped with 
the drift term on the l.h.s. of the kinetic equation. As for the identical spin rotation effects, 
they average to a spin molecular field which is the equivalent of the molecular field considered 
by Silin 3 and Leggett j^; the ISRE is already present to zero order gradient expansion, as 
opposed to the scalar retardation terms jl3[ I15j . which explains why the dispersion relation 
of spin waves is different from the usual sound wave dispersion. A detailed study of forward 
scattering in a binary collision can therefore lead to a microscopic understanding of the origin 
of the Landau mean field (lateral scattering introduces collisional damping) . From a practical 
point of view, it remains true that the mean field treatment of the drift term, necessarily 
associated with the use of a pseudo-potential, is often more compact and elegant than a de- 
tailed study of the collisional properties in the forward direction; it is nevertheless interesting 
to check precisely to what extent they are equivalent in various situations. 

Our purpose in this article is twofold. First, for particles with spins, the equivalence in 
question has been verified in the literature only to zero order in the gradient expansion, in 
other words only for the local terms appearing in the kinetic equation. Here we wish to 
derive the kinetic equation to first (non-local) order in the gradient expansion. Moreover, we 
will consider the situation where the scattering may depend on the relative spin orientations 
(even if this dependence was relatively weak in the experiments of the JILA group Ej). 
We will therefore examine in detail the properties of collisional interference in the forward 
direction in this more general case; the details of the calculations are given in two appendices 
(gradient expansion of the interaction term within the S'-matrix collision Ansatz). Second, we 
also wish to study situation where the gas is far from local equilibrium. In fact, most of the 
work on spin waves in the literature deals with hydrodynamic situations where, at every point 
of space, the gas is close to local equilibrium, so that a simple hydrodynamic expression of 
the spin current can be used; of course, this reduces the number of variables, but at the price 
of an approximation which is not necessarily justified. In fact, in the experiments in question, 
the gas sometimes evolves very far from an hydrodynamic regime, which is not so surprising 
since the spin waves have a large amplitude and lead to an almost complete segregation of 
the two spin species. For a study of spin waves in the collisionless regime (mean-free-path 
large compared to characteristic lengths), but still close to equilibrium (linear regime), see 
Ref. [lEl. 

2 Kinetic equation 

We study an ensemble of identical atoms with two internal levels, obeying either Bose or Fermi 
statistics. As mentioned in the introduction we can, without any loss of generality, assimilate 
these two levels to two spin levels, whatever their real physical origin is. For instance, hyperfine 
atomic structure states, even in a situation of intermediate magnetic decoupling, are possible; 
the only important thing is that the quantum states describing the internal variables should 
be orthogonal. We now derive a kinetic equation for a dilute gas of such atoms. 

2.1 Collisional methods; Yvon-Snider and Lhuillier-Laloe equations 

In classical statistical mechanics, Boltzmann's intuitive method for deriving a kinetic equation 
is based on the study of individual collisions; one adds the effect of all possible binary collisions 
taking place at each point of space on the single particle distribution in phase-space /i(r, p). 
Another point of view starts from the infinite hierarchy of BBGKY equations and closes its 
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first equation, which relates fi to the two particle distribution function /2, by studying the 
evolution of /2 during a binary collision and expressing it as some functional of /i (see for 
instance Ref. [T7|[TH1). 

The same lines can be followed in quantum mechanics, but /i(r,p) has to be replaced 
by the reduced single-particle density operator pi. Since we assume here that the particles 
have two internal states, which we treat as spin states, the single-particle density operator pi 
acts in the product space of orbital and spin variables of one particle. The operatorial kinetic 
equation has the form: 

^n, 1 r . 1 rl 

PI (1) 



dpi 
dt 

where: 
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//i = i + y^^* (2) 



is the single-particle Hamiltonian, p the momentum of the particle, m its mass, 1 the unity 
operator in spin space and y^^^ the operator describing the external forces acting on the 
particle, which may be spin dependent. The r.h.s. of eq. takes into account the effects of 
binary interactions between the particles; it contains the result of the collisional approximation 
mentioned above, and will be discussed in more detail below. The equation is only valid for a 
dilute gas (n"^/^ ^ a, where n is the number density of the gas and a the scattering length) 
and on time/length scales much greater than the duration of a collision/scattering length. 

The next step is to introduce the Wigner transform ^1120] pw i^,p) of pi with respect 
to orbital variables: 

pw{r,p)^{27Th)-^ I dVe^P-'/^ (r-^|pi|r+^) (3) 

It is a classical function of position and momentum, but still an operator - or a 2 x 2 matrix - in 
spin space. In the following, we will treat the orbital degrees of freedom semi-classically, but 
not the spin degrees of freedom. For the sake of simplicity, from now on we will drop the index 
W, assuming that the dependence on r and p is sufficient to signal a Wigner distribution. 
One then performs a Wigner transform of equation which involves taking the Wigner 
transform of products of operators. This is possible by using Groenewold's formula j2H I2()j. 
which provides the result as a infinite series of gradient expansion: the first term is merely 
the product of the Wigner transforms of the two operators, followed by product of gradients 
with respect to r and p, followed by higher order gradients, etc. The Wigner transform of 
the l.h.s. of equation then gives: 

dtp{r,p) + P • Vr/5(r,p) + l[p(r,p),F-*(r)]_ 
m in 

-i[Vp/>(r,p),.Vry^^*(r)]+ + ... (4) 

The first two terms are exact, whereas the next two are only the first orders terms in a 
gradient expansion involving the external potential. The commutator is obtained to zero 
order of the gradient expansion while the anticommutator occurs only to first order, followed 
by higher order terms symbolized by the dots. The anticommutator corresponds to the 
classical force term (describing the effect of a trap, for example), which we wish to retain in 
our calculations; as a consequence, we have to include up to first order gradient expansion in 
the Wigner transform of the r.h.s. of equation . 
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The non-trivial part in the derivation of a kinetic equation begins when one gives a precise 
expression to the formal r.h.s. of equation This is the subject of many studies in the 
literature and even the subject of books; here, we limit ourselves to a simplified discussion of 
two different approaches which are appropriate for the study of spin waves. We consider only 
collisions taking place in elastic channels, where the number of atoms in each internal spin 
state is conserved, and begin with some simple considerations on elementary collision theory. 
For each channel, at very low temperatures, collisions occur in the extreme quantum regime 
where the typical s-wave scattering length a is small compared to the de Broglie thermal 
wavelength Ay. In this case, collisions are well described by the isotropic s-wave, and the T 
matrix can be approximated by: 



(terms of order should be retained in order to satisfy explicitly the unitarity of the S- 
matrix in the optical theorem); in this equation, m is the mass of the particles, /fc(^) the 
scattering amplitude which, for s-wave scattering, is independent of the polar angle 6, and k 
the collision wave vector. 

A first obvious remark is that the intensity of the spherical scattered wave is proportional 
to the square of a, while interference effects in the forward direction between the incoming 
plane wave and the spherical scattered wave are proportional to a itself. Generally speaking, 
the relative magnitude of forward and lateral scattering effects is of order Xt/cl, which is 
large at low energies. This is why, here, we concentrate on forward scattering effects and 
their relation to mean field corrections to the drift terms. Another remark is that, because 
waves interfere only if they are in phase, phases and i factors are important here. It turns out 
that, in quantum mechanics as well as in optics |22[ I23j . the summation over many scatterers 
in the forward direction introduces an additional i factor (which is actually the origin of 
the i factor in the optical theorem). Hence, in every scattering channel, equation @ shows 
that no first order a interference effect occurs in the forward direction, the first contribution 
arising from the second term in ika^; to first order in a, only the phase of the forward wave 
is changed, not its intensity. But, if several scattering channels are open, the situation is 
different, because these different phase shifts in the forward scattered waves can introduce 
spin rotation. Actually, this is precisely the origin of the first order identical spin rotation 
effect in the forward direction |13| i6j . 

2.1.1 Spin-independent interactions 

We now briefiy summarize two approaches that have been developed in the context of par- 
ticles with spin: the Yvon-Snider approximation, which emphasizes more the relation to the 
BBGKY hierarchy, and closes the infinite set of equations with an approximation suitable for 
a binary collision; the Lhuillier-Laloe (or S matrix) approach, which is closer to the initial 
approach of Boltzmann and makes use of the S collision matrix. We begin with the simplest 
case, when the interactions are independent of spin. 

Yvon-Snider An interesting method is that originally proposed by Yvon |23], then in- 
dependently by Snider 2^; it provides a result which, in the literature, is often called the 
Waldmann-Snider equation; the main idea is to express the collision term as the trace over a 
collision partner (particle 2) of a commutator containing the binary interaction potential V12 
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and a product of single particle operators modified by the unitary transformation associated 
with the M0ller collision operator $7 = ri^+^i 



di 



coll 



pi(l) = (ih)-^ Tr2 \ \Vi2 , n /5i(l)/>i(2) 17t 



(6) 



This equation is valid for dilute gases of distinguishable particles; for non-degenerate gases 
(nAj. <C 1, where At is the thermal wavelength of the particles) of bosons or fermions, it 
generalizes into: 



d_ 

dt 



pi(l) = (ih)-' Tr2 



1 + ePe. 



Vi2 , npiii)pi{2) 



1 + ePe: 



coll I 

where Pgx. is the exchange operator between particles 1 and 2, and where: 

e = 1 for bosons, e = — 1 for fermions 



V2 



(7) 



(8) 



Wigner transforms can then be applied to these expressions to complete the derivation of 
the closed kinetic equation for p{r,p). For instance, expanding (jUJ in a Groenewold gradient 
expansion of the products provides, to zero order in gradients, the straightforward Boltzmann 
expression for the collision term: 



di 



/o(ri,Pi 



coll 



dh^ 

m 



crxik) p(ri,pi) /(ri,p2) 



- d k' ak{9) p(ri,p'i) /(ri,p2) 



where /(r, p) is the spin trace of p(r, p): 

/(r,p) =Trs{/5(r,p)} 



(9) 



(10) 



and where k is the unit vector in the direction of k and 9 the angle between k and k . Here 
we use the usual notation in Boltzmann theory: 



q = 2?ik P2 = Pi - q 

= p - (q/2) + hkk' P2 = P - (q/2) - hkk' 

as well as the definition of the differential and total cross sections: 

ak{e) = {ATT^m'/ti') |r(k,k')|' ; cjTik) = fd^k' akie) 



(11) 



(12) 



Higher order terms will introduce scalar mean field terms, which we do not write here explicitly 
for concision; they can be found in the appendix of ref. p[4j. 

Now, if we introduce spin and statistics, we have to use instead of To zero order 
in gradients, this introduces the additional exchange terms to the right hand side of Q: 



-| j d^Q ^ {irf^d. (k) [p(ri, pi), ^(ri, p2)]_ 
- J d-'k' irl^ie) [/5(ri,p;), p(ri,p^)]_ + 
/ d-'k' al-iO) ([p(ri, pi), p(ri, p2)]+ - [/i(ri, p'^), p(ri, p^)]^) } 



(13) 
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where the "generalized cross sections" are defined by: 

Tf^Ak) = {-87T'm/h''k) Re{r(-k,k)} (14) 

and: 

ar-(0) - «rf'-(0) = (47rW//i^) r(-k,k') r*(k,k') (15) 

In (|13() , the two terms in the first two lines correspond respectively to the ISRE in the forward 
and lateral direction; the two terms in the third line correspond to exchange effects changing 
the values of the total and lateral cross sections, as discussed in more detail in For 
simplicity, we do not write the first order gradient terms, but more details on the limit of 
validity of the Yvon-Snider equation and of the various terms which it contains can be found 
in |26j . jl5j or jl4j (in particular the Appendix of this last reference). Inserting the low energy 
limit is trivial; we write the explicit results in the next section. 



Lhuillier-Laloe Another point of view, more directly in the spirit of the initial Ansatz of 
Boltzmann, was employed by Lhuillier and Laloe (LL) it is based on the use of the S 
collision matrix. When collisions are treated only as "closed processes" (one ignores particles 
"in the middle of a collision"), it is indeed possible to use the S matrix to relate exactly the 
single particle operators after and before collision. One gets the following expression for the 
single particle density operator after collision as a function of the same operator pi before 
collision: 

P'i(l) = Tr, S PI (2) St i^^l (16) 

(the trace acts over orbital and spin variables of the collision partner 2) where the S matrix 
is related to the collision T matrix by: 

7TTYI 

(1 : k/; 2 : -k/ I 5 I 1 : ki; 2 : -k^) = <5(k/ - k^) - ij^ 6{kf - h) T^, (17) 

{hk is the relative momentum of the two particles) . One then approximates the rate of change 
of pi by the variation {p[ — pi)/ At during a short time interval ^ At, which is however much 
longer than the duration of a collision, to obtain: 



It 



1 _ \l + ePe: 



pi(l) — I -^-jJ^ [S p,{l) p,{2) - p,{l) pi(2)J -^-jJ^ \ (18) 



1 + ePe: 



One should note that the S matrix Ansatz of equation (|18j) is only valid for a non-degenerate 
gas where: 

l^^-^Pi(l) Pi(2)^^^| = + ehilf - Pi{l) (19) 

As above, the kinetic equation is obtained from equation (|18j) by performing a Wigner 
transform followed by a gradient expansion. In Appendix I we first calculate the zero order 
gradient expansion, which is given by equation (|7H|) ; it is easy to see that it coincides exactly 

^We will see that At is a time which emerges naturally from the calculation through the occurrence of 
square of delta functions of the energy (see Appendix I) . 
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with the sum of © and which shows that the Yvon-Snider and 5 matrix approximations 
are exactly equivalent to this order. The first order gradient terms are also given in the same 
Appendix; here, we only give their low energy limit by using equation @, which provides the 
following kinetic equation: 



5tp(ri, pi) + ^ • Vr/5(ri, pi) + 4 kn, Pi), V'^^'iri) +egh{ri) 
m in I 



I coll 



[p] 



(20) 



- 2 [Vp;5(ri, pi), •Vr(y^"*(ri) + gn{ri)i + egh{r,)) 

where g = Anh'^a/m and n(ri) is the local density operator integrated over velocities: 

n(r) = J cfppir,p) (21) 

In this equation, as usual, forward scattering terms linear in the coupling constant g oc a 
have been included in the l.h.s. The local term of the collision integral, which regroups terms 
proportional to a^, is given by: 



^ colli 



dq— / dk'a /5(ri, pi)/(ri, P2) - /5(ri, p'i)/(ri, ps) 



m 



+ 2 [p(ri,Pi),p(ri,P2)]_ 



[/5(ri,p'i),/5(ri,p^)] 



(22) 



where p2, p'l and P2 are defined in (jllj) . In addition to the anticommutator in the l.h.s. 
of (|2()j) . obtained to first order in the gradient expansion, we recover exactly the low-energy 
limit of the LL equation pH] with, in the l.h.s., the ISRE in the forward direction (the ISRE 
in the lateral directions disappears in the low energy limit); the second line of (|22|) provides 
the statistical terms introduced by particle indistinguishability into the total and differential 
collision cross section. 



2.1.2 Spin-dependent interactions 



We now study the gas with spin-dependent interactions and, for a moment, consider the 
particles as distinguishable. In the limit of slow elastic collisions, we need only consider four 
s-wave scattering lengths: an, 022, a^^2 '^12 • '^^^ ^^^^ describe collisions processes 
between particles that are in the same internal state; when they are in othogonal internal 
states, two different processes occur for distinguishable particles: a direct collision process 
without energy transfer described by , and a transfer collision process described by " 
similar processes occur in the theory of spin exchange collisions |27j . The T matrix is now a 
4x4 matrix in spin space; to lowest order in the scattering lengths, its matrix elements with 
distinguishable particles are given by: 

/ an (A;) 




(l:a;2:/3|rfc| 1:7;2:5) 



27r2 



m 



\ 















'12 



Ad) 

n2 



(23) 



(k) J 



where a, /3, 7 and S are the internal state indexes taking values {1, 2} and where: 

aa/sik) = a^p [1 - ikaajs] 
Equation (|23|) replaces Q when the interactions depend on spin. 



(24) 
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Kinetic equation The derivation of the kinetic equation then proceeds along the same 
hnes as in the previous section: we use ()16p again to symmetrize the density operator and 
take quantum statistics into account, but now insert (|23|) instead of (|17|) in it, calculate 
the partial trace, and obtain an interaction term for the single particle density operator; 
we then Wigner transform the result and use a gradient expansion limited to first order. 
This calculation is similar to that of Appendix I for spin-independent interactions, but of 
course more complicated; it provides more general expressions for the mean field terms and 
the collision integrals. For conciseness, we limit ourselves to the mean field terms, which will 
actually be sufficient for our calculations below. In this particular case, the terms can actually 
be obtained more conveniently by an operatorial method, different from that of Appendix I, 
as shown in Appendix II. We only reproduce the final kinetic equation: 

dtp{r, P) + - • Vr/5(r, p) + hp{r, p), U (r)]_ 
m in 

- i[Vpp(r,p), •Vrf/(r)]+ = I,ou[p] (25) 
where the effective single-particle potential U is: 

fj^yext^f (1 +e)5'22n2 +512^1 £912^21 \ ^26) 

V egi2ni2 (1 + e)giini + 512^2 / 

with gap = A-Kh^aafj/m, where 012 = 0^2^ +ea^2 ; it turns out that only this combination of the 
direct and transfer coupling constants enters the calculation (for more details, see Appendix 
II). 

This kinetic equation resembles the Landau-Silin equation j31|2H] for a normal Fermi liquid 
in a magnetic field. The second term in the l.h.s. of (|25() is the usual free drift term; the 
anticommutator is a force term including both the effect of an external potential and of the 
mean field; the commutator is a spin precession term containing effective magnetic fields from 
different origins (differential Zeeman effect, spin mean field, etc.). 

Evolution of the density and spin density in phase space The structure of the 
distribution function, which is still a matrix in the space of internal variables, is more explicit 
if we decompose p and tJ in the basis formed by Pauli matrices a and the unit 2x2 matrix: 

p(r, P, i) = ^ (/(r, P, t)i + M(r, p, t) ■ a) 
U{v,t) = Uo{v,t)i + \]{r,t)-a (27) 

When written in terms of the phase-space density / and spin density M, equation H25() 
becomes: 

dtf + ^- Vr/ - VrC/o • Vp/ - V,U • VpM = M] 

OT T 

dM + ^- VrM - Vr?7o • VpM - VrU • Vp/ - — X M = I^^l} [/, M] (28) 

In these equations, Uq plays the role of the overall trapping potential and U is the local 
effective magnetic field around which the spins rotate. The effective trapping potential Uq is 
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given by: 



jrext I T/ext 



1 + e 
2 



(522 +511) +512 



n 
2 



1 + e rriii 
+ ^(522-911)^ (29) 

and the effective magnetic field is: 

U(r, t) = ^^^^'^^ c|| + c ^i2m(x,t) 

;il^(r, t) = 1^2-* - + ^ [(<722 - gii)n + (^/22 + 511 - 25i2)m||] 
where the local density and spin polarization density are defined as: 



(30) 



n = 71-2 + ni ; my = n2 — 7^l 
m-x.i = "-21 + ni2 ; m_L,2 = ^("-21 - "-12) (31) 

The basis in spin space is denoted by {e_L^i; e±^2', e\\}- The two first vectors define the 
transverse plane, the last defines the longitudinal direction. 

The collision integral in (|25|) . i.e. the equivalent of eq. (|22|) when the interactions are 
spin-dependent, is not written explicitly; to simplify the calculations, and because we will 
mostly study experimental situations where lateral collisions do not play a dominant role, we 
limit ourselves to a simple relaxation time approximation: 

(/) /(r,p,t)-r'^(r,p,t) ^(M)„ M(r,p,t)-M^'?(r,p,t) 

^coll — ^ ' ^coll ~ ^ V-J^j 

where /^"^ (resp. M'^'') is the local equilibrium (spin) density in phase-space, and r is a 
relaxation time. Actually, one could allow for different relaxation times for /, My and M_l. 
Nevertheless, in the case of a non-degenerate gas with spin-independent interactions, it is 
known (Ref. [211 and references therein) that these relaxation times are equal; below we 
consider a non-degenerate gas of ^''Rb, for which the three scattering lengths an, 022 and 012 
are very close PP, so that in practice we will we only need a single relaxation time r given 
approximately by: 



1 / m 



(33) 



47raf2 n y ksT 

Finally, the kinetic equations we will use below to discuss spin waves are: 

dtf + -- Vr/ - Vrf/o • Vp/ - VrU • VpM ~ -l^IH 

m T 

p 2U M - M'^* 

dtM + — • VrM - VrUo ■ VpM - VrU • Vp/ - — X M ~ (34) 

m n T 

The main differences with the LL kinetic equation derived in Ref. |13j are that the full mean 
field is included (not just the spin mean field), and that the equation is obtained for general 
spin-dependent interactions. Nevertheless, our result is also less general, since the collisions 
are assumed to occur only in the s-wave channel and that the collision integral is treated at 
the relaxation time approximation level (if necessary, these two restrictions could be lifted 
without any special difficulty). 
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2.2 Mean field 



Another way to derive the modified drift term is to directly use a mean field approximation 
within field theory. The following second quantized Hamiltonian density is used to describe 
a trapped ultra-cold gas of bosons or fermions: 

H =Hl + Hint 

= E (|^VVi(r)-VV.(r)+V;-*(r)V^t(r)v,„(r)') 
a=l,2 ^ 

+ ^ j^^^/^V'iw^w^/^wv'aW (35) 

The annihilation and creation field operators of a particle in the internal state a are '0o(r) 
and ipa{r). They have bosonic or fermionic equal-time commutation or anticommutation 
relations. The interactions only depend on the three coupling constants gu, gu and g22', the 
real potential ^-matrix elements have been replaced by the appropriate T-matrix elements at 
low-energy, according to an usual procedure in the study of cold gases (Fermi pseudo-potential 
method or V ^ T renormalization procedure in Ref. j2Sj). 

We wish to obtain the time evolution of the single-particle density matrix pap{r,T',t) = 
{ijja{Y' , t)%l)p{r, t)), where (. . .) denotes the expectation value with respect to the time-dependent 
A'"-body density operator of the system. To this end we write down the Heisenberg equation 
of motion for '0a(r', f)^/)^(r, t) using the Hamiltonian H35|) . The non-interacting part of the 
Hamiltonian density TCi gives the single-particle Liouville-von Neumann equation: 

ihdtpai3{r, r', t) = -— (Vr - Vr/) Paf^{r, r' , t) 

+ (yr*(r)-^r*(r'))p./3(r,r',i) (36) 

The interacting part Hint couples the single-particle density matrix with the expectation 
value of the product of four field operators (two-particle density matrix). If we are only 
interested in time scales large compared to the duration of a collision, we can decompose the 
averages of such products into averages of products of two operators ^ . In the context of the 
equation of motion for the single-particle density matrix, this approximation is often called 
the random phase approximation (RPA) j31j . The interacting part of the equation of motion 
then becomes: 

E C^(r,t)p5/5(r,r',t) - ^ p^s{r,r' ,t)V,^/ {r' ,t) (37) 

(5=1,2 (5=1,2 

where the mean field potential 2x2 matrix is given by: 

^r/i^) = E 5a7(V'^(r)V'7(r)) + ego.p{^l,i{v)i,p{v)) (38) 

7=1,2 

which contains a direct term and an exchange mean field term, proportional to e. Within 
this approximation, each particle evolves under the infiuence of an effective single-particle 

■^This result can be obtained for instance from the Wick theorem. It is valid because the gas is not Bose 
condensed; otherwise, the mean field could not be obtained in that way, due to the absence of the exchange 
term. 
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potential U{r,t), which is the sum of the external potential y^^'*(r) and the time-dependent 
mean field V"^-^ (r> 0- 

In order to compare this result with that of the preceding section, we note that ^^^(r, t) = 
{ipa{r, t)?/'Q,(r, t)) is the average density of atoms in internal state a = 1,2 and the off-diagonal 
elements ni2(r, t) = n2i{r,t) = {^p\{r,t)'^lJ2{r,t)) are the coherences between the internal 
states, so that the effective potential U reproduces exactly the one obtained in the preceding 
section, equation ()26() . We therefore arrive at the following operatorial kinetic equation: 

= (39) 

which is the quantum equivalent of the Vlasov equation (see Ref. jl8| for example); in this 
equation. Hi is the single particle hamiltonian defined in ((21) and V'^-^ is defined by (|38|) . 
Obviously, if we Wigner transform this operatorial equation, we will obtain a kinetic equation 
that is equivalent to (|25|) . but without the collision integral. The effect of lateral collisions can 
then be treated phenomenologically by adding by hand a term corresponding to a relaxation 
time approximation, as usual then equations (|3i|l are recovered and the mean field 
treatment of interactions becomes exactly equivalent to these simplified equations. 

2.3 Discussion 

We have used three different methods to derive a kinetic equation for a gas of atoms with 
two internal levels: the S matrix (or LL) Ansatz, the Yvon-Snider Ansatz (YS) and the mean 
field approximation; we now summarize a comparison between the results. 

To zero order in the gradients, LL and YS give exactly the same equation, containing a 
Boltzmann-type collision integral and a spin mean field. As for the mean field approximation, 
it provides only the spin mean field, obtained in the limit of low-energy collisions because of 
the use of the pseudopotential; this limit is in full agreement with the results of the two 
other methods. We note in passing that this spin mean field occurs for a non-condensed 
Bose gas, but not in a Bose-Einstein condensate (BEC) at very low temperature, because 
the exchange term is absent when a single quantum state is involved. In this case, there is 
however a different mechanism for creating spin waves, the so-called "quantum torque" of 
purely kinetic origin [32] ; it arises naturally because, in quantum mechanics, any gradient 
of the phase corresponds to a kinetic energy for each component which, in turn, affects the 
evolution of the relative phase and therefore of the transverse spin orientation. 

To first order in the gradients and for forward scattering, LL, YS and the mean field 
approximation all give the same scalar mean field term (zero sound), again in the limit of 
low-energy collisions for the third method. We have checked that the complete expression of 
the LL and YS first order gradient terms coincide exactly for forward scattering, but not for 
lateral scattering; whether or not the two methods are completely equivalent to this order is 
left as an open question; for more details see |33] . 

In a non-degenerate dilute gas, the spin mean field generally dominates over the scalar 
mean field, because it appears to lower order in the gradient expansion; it is therefore not 
surprising that zero sound type collective modes should not propagate as easily as spin waves 
in dilute non-degenerate gases. 

To what extent collisional methods and mean field approximations are equivalent is not a 
trivial question. Generally speaking, mean field theory is assumed to be valid for systems with 
either long range interactions (Vlasov's equation for a plasma JHl) or in strongly interacting 
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degenerate systems described in terms of quasi-particles (Landau's equation for a Fermi liquid 
for a dilute gas, neither of these conditions is met. Our conclusion is nevertheless 
that the equivalence is perfect to first order in the scattering length a but does not hold to 
higher orders in a. We finally remark that, when describing a non-degenerate dilute atomic 
gas in terms of quasi-particles (see the work of Bashkin ^B] , for example) , a quasi-particle is 
just a particle whose kinetic energy /2m is shifted by a local mean field term, so that the 
kinetic drift term of the quasi-particle transport equation is left unchanged. 

3 Spin oscillations and internal conversion 

Our purpose now is to use the kinetic equation obtained in the preceding section to discuss spin 
waves in ultra-cold trapped atomic gases. We first briefly review a recent experiment realized 
at JILA 1 and the theoretical work stimulated by it [HI, 35, 36||SI1ISH]- We then discuss the 
non-hydrodynamic character of the observed spin waves, which allows a comparison between 
the different theoretical treatments. Finally, we describe the decay of the spin waves and 
discuss the relevance of decoherence. 

3.1 The JILA experiment 

In a recent experiment performed at JILA , Lewandowsky et al have studied bosonic ^^Rb 
atoms with two hyperfine states of interest (denoted by 1 and 2), confined in an axially 
symmetric magnetic trap elongated along the Ox direction. Their experiment is described in 
PP; the temperature T of the gas in this experiment is about twice the critical temperature 
for Bose-Einstein condensation, so that the gas is not strongly degenerate and can be treated 
reasonably well as a Boltzmann gas; however, since the de Broglie thermal wavelength is 
much larger than the average scattering length, collisions occur in the full quantum regime. 
Initially, the gas is at equilibrium with only state 1 populated. A 7r/2 radio frequency pulse 
is then applied, which suddenly puts all the atoms into the same coherent superposition of 
states 1 and 2. The subsequent evolution of the system along the axial direction is then 
monitored by measuring optically the local densities ni and 77-2 of atoms in each internal 
state. Experimentally, one observes that the system "segregates" : after about 100 ms, atoms 
in state 1 are mostly found away from the center of the trap, while atoms in state 2 move 
towards the center, the total density n = ni + n2 remaining practically unchanged. After 
about 200 ms, the local densities of each species return to equilibrium. 

In Ref. 1^, this phenomenon was explained in terms of an internal conversion resulting 
from the identical spin rotation effect (ISRE), which is at the origin of the spin mean field; 
similar considerations were almost simultaneously provided by two other groups |34l I36j . 
Qualitatively, this spin oscillation can be understood as follows: 

(i) the field gradient creates an inhomogeneous spin precession, so that the gas develops 
a gradient of transverse spin polarization 

(ii) the thermal motion of the atoms then creates correlations between transverse spin and 
velocity; therefore, a particle moving with a given velocity at point x gets a spin polarization 
which is not parallel to the average local polarization 

(iii) the ISRE then makes its spin polarization leave the transverse plane and get a longi- 
tudinal component, in a direction which depends on the sign of the velocity 

(iv) this is equivalent to a velocity dependent internal conversion which, at some later 
time, results in an efficient spatial separation of atoms in the two internal states. 
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To study this phenomenon quantitatively, we derive from (|34() an effective one-dimensional 
kinetic equation in terms of the local density in phase space f{x,p,t) and spin density 
M(x,p, f), by assuming that radial equilibrium is quickly established (the radial trap fre- 
quency a;rad/2vr is much larger than the axial trap frequency LOax/'^T^)- Integrating 1)34(1 over 
radial coordinates y^z and momenta Py,Pz provides the following equations: 

dtf + -dj- dxUo dpf - dxV ■ dpM ~ -(/ - D/t (40) 
m 



^ d^m - dxUo dpM - a,u dpf-^xM 

m h 

~-(M-M'^9)/r (41) 



(every quantity depending initially on three-dimensional coordinates r and momenta p has 
been integrated over radial coordinates and momenta, so that it now depends only on the 
coordinate x and the momentum p = Px)- In the process of radial averaging, the coupling 
constants ga/s are renormalized by a factor 1/2, as discussed by Levitov JSl!- The three- 
dimensional effective trapping potential ()29|) therefore becomes: 



\rext _|_ j/ext 



1 + e 
2 



(522 +511) +512 



n 
4 



1 + e mil 
+ ^(522-911)^ (42) 

and the three-dimensional effective magnetic field ()3()|) is changed into: 

U(x,t) = ^^e||+e^l^^^ (43) 



with: 



1 + e 



2 



n , mil 



hn{x, t) = - + {922 - 5ii)t7 + (522 + 511 - 2512) 



2 " " ^ 2 

where ey is the unit vector in the longitudinal direction in spin space. 

A few simplifying assumptions are appropriate for the experimental conditions of Ref. 
The confining energy (V2^^* -|- Vi^'')/2 is of order /cbT ~ 13 kHzx/i which is much larger 
than the mean field interaction energy gn{0) ~ 140 Hzx/i. This allows us to keep only the 
confining energy of the harmonic trap in the effective trapping potential Uq: 

^„(,,.S3^)+S!!(£).^„,.i.. ,44) 

The differential trapping energy V2^^ — Vf^^ ~ 10 Hzx/i is even weaker than the mean field 
interaction energy. This ensures that the force terms dxV-dp{f or M) appearing in equations 
()40I41|) are negligible. On the contrary, when the local effective magnetic field U does not 
appear under a spatial gradient, as in the term 2U x M/^, it can not be neglected. This 
effective magnetic field U is made of two terms (see equation (|43() ): one is an effective external 
magnetic field Ml/ 2] the other is an exchange magnetic field or spin mean field (7i2ni/4 which 
results from the ISRE. The effective external magnetic field is the sum of the contributions 
of a differential Zeeman and of a differential mean field: 

hn{x) = h n{x)eu ^ [y-*(x) - Vr\x) + (522 - giiMx)/2]eu (45) 
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where, following ref. pQ, we have assumed that 2gi2 — 511 + 522 for simplicity. The average 
value over the sample of the effective external magnetic field can be removed by going to a 
uniformly rotating frame (Larmor frame). For the numerical simulations, we need to know 
the X dependence of this effective magnetic field. As the profile n{x) will be shown to be a 
Gaussian (as for a non-interacting non-degenerate gas) which does not vary in time, we model 
the experimentally measured effective external magnetic field by: 

n{x) = -Sn exp -(ma;2^xV2A;BT)e|| (46) 

where the parameter 6Vt is the variation of Vt{x) between the center and the edge of the cloud. 
Taking into account all approximations mentioned above, we finally arrive at equations: 

dtf + mcolxdj,/ -(/ - D/r (47) 

dtM + — d^M - mujl^x 9pM - iVt + ^i^) x M 
m In 

~-(M-M^«)/t (48) 

The initial equilibrium Maxwell-Boltzmann distribution f{x,p) solves the kinetic equation 
(|47|). so that the dynamics after the 7r/2 pulse can be expressed in terms of M only. In 
addition, the density n(x) is obtained by integrating the equilibrium Maxwell-Boltzmann 
distribution over momentum p, which shows that it is time independent; from now on, we 
will only consider equation (|4(S|) . The spin distribution immediately after the 7r/2 pulse is 
assumed to be: 

M{x,p,t = 0) = f{x,p)e±^i (49) 

Equation l|48|) can be solved numerically with values of the parameters taken from Ref. 
^: the trap frequencies are iVax/'^'^ = 7 Hz and uj^ad/'^'^ = 230 Hz; the average time between 
collisions is r ~ 10 ms; the temperature is T ~ 0.6 ^K; the density at the center of the 
trap is n(0) = 1.8 x 10^^ cm~^; 6Q/27T is typically ~ 12 Hz, and the scattering lengths are 
Oil = 100. 9ao, a22 = 95.6ao and 012 = 98.2ao where oq is the Bohr radius. The solution at 
the center of trap is plotted in Figure it shows very good agreement with the experimental 
observations, without any adjustable parameter. 



3.2 Large, non- hydro dynamic, spin waves 

Ordinary spin waves in non-degenerate gases usually occur either in the hydrodynamic regime 
or / and in a regime of small amplitudes |H1 El ^1 ■ One can then assume a small departure of 
the spin distribution (from either local or global equilibrium); this is similar to the first order 
Chapman-Enskog gradient expansion, which leads to the Navier-Stokes equations (see for 
instance pUjl ) . The result is the well known Leggett equations |1] , transposed from degenerate 
Fermi liquids to dilute non-degenerate gases P: 

(9tm + (9J = n X m 
9tj-(n + ^)xj + %fa.m + a;L^m -1 ^"^^ 

where the spin current along Ox is defined by: 

jix,t) = [dp ^M{x,p,t) (51) 
J m 
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Figure 1: Time evolution of the spin polarization m at the center of the trap (x = 0) when 
517/27r = 12 Hz. 

The momentum dependence has disappeared from these equations, which contain only po- 
sition variables. These equations can describe spin waves, not only in the hydrodynamic 
regime, but also in the collisionless regime as, for example, in liquid '^He (Silin spin waves 
|39j ) or in H| gas |4f)) . For a discussion of the hydrodynamic-like description of spin waves in 
the collisionless regime of non-degenerate dilute gases, see Refs. [XHIIIT] . 

The spin waves observed at JILA were of large amplitude and occurred in the intermediate 
regime between hydrodynamic and collisionless, where u^axT ~ 1, so that there is no a priori 
reason to believe in the validity of the Leggett equations; we therefore need to resort to a 
numerical solution of the kinetic equation in terms of both position and momentum variables. 
Figure [21 shows the results of this calculation, and the spin density distribution (longitudinal 
and transverse) in p-space for different times. For comparison, the local equilibrium spin 
distribution is also plotted (dashed line); we see that the spin distribution in p-space can 
indeed get very distorted, which illustrates a strongly non- hydrodynamic situation. For very 
short times, t < 4Q ms, the transverse and longitudinal spin distributions are still close to 
local equilibrium {M^^ = 0), so that the Leggett equations are valid. For t > 40 ms, this is no 
longer true, so that one would lose a significant part of the physics by not retaining the full 
momentum dependence of the distribution. In other words, studying the two first moments 
of the distribution (m and j) is not equivalent to studying the full distribution M(p). What 
actually takes place is a phenomenon analogous to a density "shock wave" rather than a 
"sound wave". Eventually, for t > 160 ms (end of the "segregation"), since equilibrium is 
almost reached, the Leggett equations become valid again. 

The initial part of the internal conversion (or state separation) can be treated analytically 
in several approximations, allowing for comparison between different approaches such as those 
of Refs. In Ref. we solved the kinetic equation analytically, for times t <C r, 

by expanding in Taylor series in time and obtained: 



where Vt' and Vt" are the first and second spatial derivative of Vt{x). This formula predicts an 
initial quartic in t behavior, which correctly reproduces the numerical results; each power of 
t corresponds to one of the four physical processes (i) to (iv) described in § 13.11 




(52) 
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Figure 2: Longitudinal spin distribution M||(x = XT,p,t) and transverse spin distribution 
-^±,1(2^ = XT,P,t) at position x = x-p for t = 10; 40; 85; 160 ms, where xt = y^tBT/Tnuj^ 
and pt = \/mkBT. The local equilibrium spin distribution is plotted with a dashed line. 
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Figure 3: Beginning of the internal conversion at the center of the trap. The numerical 
simulation is plotted with circles, the result with a full line and the result with a dashed 
line. 



Another approach to calculate the short time behavior is that used by Williams et al [HE] • 
They use the Leggett equations which are valid for t <C r because the spin distribution is at 
equilibrium at i = 0. Solving the Leggett equations in the small time limit, where — j/r is 
negligible, they recover ((52). 

A third approach is to consider, following Oktel and Levitov |34j ) . the Leggett equations 
in the hydrodynamic regime. The approximation consists in assuming that the spin current 
remains close to its stationary value, which allows to neglect the term in equation (|5U() . 
Using the resulting equations, we obtain a different behavior for small times: 

n[x) In m 4! 

(the effect of the hydrodynamic approximation is to replace by tr). The and results 
are plotted in Fig. |31 we see that times smaller than a time between collisions r can not 
be dealt with in this approach. Nevertheless, after each atom has made one collision on 
average, the spin current is close to its stationary value and the hydrodynamic approximation 
is reasonably valid. Therefore, when t > r ~ 10 ms, we get: 

^ ~ _ (54) 

n[x) In m 4! 

(the parameter to includes the accumulated effect of the retardation of the spin current j for 
small time). The {t — to)^ behavior is plotted in Fig. \^ (the reason why, after 40 ms, the 
numerical solution of the kinetic equation is not well approximated by the (t — to)^ law is that 
the spin distribution gets very distorted, so that the Leggett equations are no longer valid). 



3.3 Decay of the phenomenon 

Up to this point, we have studied mostly the initial part of the phenomenon, when the 
amplitude of the spin wave grows and the two species tend to separate from each other; 
we now discuss its final part, when the spin oscillations decay and, eventually, the system 
returns to equilibrium. A first remark is that it is relatively easy to obtain an order of 
magnitude of the time at which the transition between the two regimes occurs, i.e. the time 
at which the maximum of the separation takes places. For short times, we have seen that 
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Figure 4: The numerical simulation is plotted with circles, the result with a full line and 
the {t — to)^ result with a dashed line (to ~ 2 ms here). Note the good agreement between 
the {t — to)^ law and the kinetic equation in the range r ~ 10 ms< t ^ 40 ms. 

this separation is an indirect result of a differential precession of the transverse component of 
the spin at different points along the Ox axis; the thermal motion of the atoms then creates 
correlations between the transverse orientation of the spins and the sign of the x component 
of the velocities. But, for long times, the differential rotation of the spins will be so large 
that each sign of this velocity will become associated with a widely open fan of transverse 
orientation, with almost zero average; clearly, the selective internal conversion effect will then 
also average to almost zero. Let us consider the simple situation of a homogeneous gas in 
a box of size L with a linear gradient of external effective magnetic field O'. We focus, for 
instance, on the situation at the center of the box. At time t, atoms initially at x = v xt will 
cross this point with a transverse orientation rotated by an angle of the order of: 

n' xxxt = n' vt"^ (55) 

Therefore, at times greater than: 

1 , , 

tm 7== (56) 



V 



the transverse directions of the spins average out so that the apparent segregation effect no 
longer takes place. In the case of the trapped gas, 0' ~ (50/L and v ~ uJaxL where L is the 
size of the cloud, and we obtain: 

^"^ /rr. ^^'> 

Typical values give tm ~ 100 ms in accordance with the results of the simulation. After this 
maximum, the system tends to return to equilibrium during a time period which we call the 
decay of the phenomenon. 

Figffl shows that the decay of the longitudinal and transverse spin polarizations are sig- 
nificantly different. The longitudinal spin polarization m|j(x,t) returns to equilibrium in two 
steps: it first quickly decreases from its maximum value to almost zero in a time scale of 
the order of 50 ms; it then goes to zero much more slowly, on a time scale of the order of 
200 ms. In order to understand this, we note in Fig^that, at the maximum of separation, 
the wave packets (or the clouds) associated with each internal state overlap only little; the 
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first step is then easily interpreted as the free motion of the clouds under the restoring force 
of the trap. It takes a quarter of a period T^xl^ = I'^^ax = 36 ms to go from the maximum 
of oscillation to the center of the trap. Actually, as soon as the two clouds start to overlap, 
they interact through a repulsive mean field ( 5 > 0), so that it takes a little more time than 
a quarter of a period to reach the situation where the two clouds overlap significantly ^. At 
this point, the subsequent evolution of the system can be understood as the mutual diffusion 
of two gases ^; this happens ~ 150 ms after the initial pulse and can again be described 
by the Leggett equations (we checked on the p distribution that the system is close to local 
equilibrium). The relaxation time scale is then given by the spin diffusion time tdiff L? / D 
(where L ~ 2\JkBT jraiJl^^ is the size of the cloud and D = ksTr jra is the spin diffusion 
constant). Using the same parameters as in the simulation, we obtain tdi// ~ 200 ms, in good 
agreement with the numerical result. 

We now discuss the transverse spin polarization m_|_(x,t). As seen on Fig^ it oscillates 
many times before going to zero. After ~ 150 ms, when the Leggett equations are valid again, 
we know jfi) that spin waves occur in the hydrodynamic regime. The spectrum derived from 
these equations predicts damped transverse spin waves 6^ . The damping time is of the order 
of the spin diffusion time tdi// ^ 200 ms, which gives an order of magnitude for the time 
it takes the transverse spin polarization to go from its value at the maximum of separation 
to zero. This is in agreement with the relaxation time obtained by numerically solving the 
kinetic equation (Fig^. We finally note that the frequency of the transverse spin waves is 
essentially given by the external effective magnetic field precession frequency bQ^jln = 12 Hz. 

3.4 "Ghost" wave packet and decoherence 

We now come back to the time at which the separation is maximum. In order to simplify 
the discussion, we consider two clouds of atoms, one corresponding to state 1 and located 
around x = d/2, the other corresponding to state 2 located around x = —d/2, where d is the 
distance between the center of the two clouds; we note 5x the width of each of these wave 
packets, assuming that 5x < d. In this case, elementary quantum mechanics predicts that the 
local transverse spin polarization vanishes everywhere (this is because the local spin density 
corresponds to an operator which is local in ordinary space). Nevertheless, the superposition 
of the two wave packets remains coherent so that, rigorously speaking, the system is not 
equivalent to a classical mixture of two gases, in opposition to what we have assumed in the 
previous section. 

One may wonder how this coherence translates in terms of the Wigner distribution. The 
answer is well known: the presence of the coherence is contained in a so-called "ghost" wave 
packet, which exists around the middle point between the two real wave packets (or clouds). 
For spinless particles, it "carries interferences with it" jH], and reconstructs them as soon as 
the real packets overlap again. For particles with spin, when the two separated wave packets 
have opposite spin orientations, the ghost wave packet centered at x = appears only in the 

^That the two clouds interact mainly trough the mean field and not trough (lateral) collisions is confirmed 
by the numerical solution of the kinetic equation. Indeed, we checked that removing the collision integral does 
not change the decay time scale. The role of the collision integral is merely to damp revivals of the longitudinal 
spin waves, which are observed when the collision integral is discarded. 

^Actually, an instability of the transverse spin polarization is possible when there is a strong longitudinal 
spin polarization gradient (which is the case when the two clouds overlap again) and a spin mean field. This 
is known as Castaing's instability |42| . In Ref. |43l . we checked that it does not play a role in the experiment 
done at JILA |T]. 
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non-diagonal spin matrix elements (transverse spin components of the Wigner distribution). 
In both cases, the wave packet rapidly oscillates as a function of momentum p, taking positive 
and negative values, with an oscillation period 27rh/d; the further apart the two clouds, the 
faster the oscillation. On the other hand, for a statistical mixture of the internal states, the 
ghost wave packet does not exist. Such oscillating wave packets are well-known in the context 
of macroscopic superpositions of states |44) . 

The next question is to what extent the ghost wave packet is preserved by the time evo- 
lution of the kinetic equation. Since the Wigner formalism in itself implies no approximation 
(it is strictly equivalent to the use of operators), the question arises only because our ki- 
netic equation does not provide the exact quantum evolution: it was actually obtained from 
the operatorial (quantum) equation for the single-particle density operator by truncating the 
gradient expansion of the Wigner transform. In other words, we made a semi-classical approx- 
imation. But, for the ghost wave packet, the p oscillation period is 2'Kh/d while the width 
of the wave packet is 5x, so that the parameter of the gradient expansion is d/6x, which 
may be larger than 1: therefore, our kinetic equation does not necessarily remain correct for 
this case; the semi-classical gradient expansion may smooth the oscillations of the ghost wave 
packet, and artificially introduce decoherence. Another related source of artificial decoher- 
ence^ might come from the discretization of phase-space (necessary for the numerical solution 
of the kinetic equation), which may introduce a lattice spacing in momentum space larger 
than 6p. 

Up to this point, no real physical decoherence mechanism was included in our discussion: 
we have assumed that the system is perfectly isolated. But this is not the case in practice. 
In order to obtain an estimate of the physical decoherence time, we proceed by analogy with 
a well-understood situation in quantum optics: a macroscopic superposition of two coherent 
states of the electro-magnetic field in a cavity. In that case, it is known that the coherence 
is lost as soon as one photon on average escapes the cavity and goes in the environment 
We therefore assume here that for our trapped gas, a source of decoherence comes 
from the fact that atoms can leave the trap, because of 3-body collisions for example. The 
time it takes one atom to leave the trap is l/a3f,n(0)^A^ where a^b ~ 4 x 10^^^ cm^/s is 
the 3-body recombination rate constant [ITj; ^ ~ 10^ is the total number of atoms and 
n(0) ~ 2 X 10^^ cm~^ is the density at the center of the trap. This gives ~ 0.1 ms as a 
rough upper bound to the physical decoherence time, which is still much shorter than the 
characteristic time scale of the spin state segregation. Our conclusion is therefore that, when 
the wave packets recombine, coherence is not likely to play a role in these experiments; this 
justifies our discussion in the preceding section. 

4 Conclusion 

We have derived a kinetic equation for a Boltzmann gas with two internal levels using dif- 
ferent methods. A first conclusion which emerges from this work is that, for dilute gases at 
low temperature, several significantly different theoretical approaches lead essentially to the 
same result. Here we have discussed the Yvon-Snider method, with its peculiar way to close 

^Concerning the second source of artificial decoherence, we think that it is not a problem in our simulations. 
We checked that, when removing the collision integral so that the kinetic equation becomes reversible, the 
transverse spin polarization disappears when the two wave packets do not overlap but reappears when they 
mix again. 
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the BBGKY hierarchy, the S matrix method which is very close to the spirit of the original 
Boltzmann equation, as well as the popular mean field method. In the latter case, the equiv- 
alence actually holds only in the low temperature limit, when forward scattering dominates 
over lateral scattering by a factor of the order of At/o- Whether or not it is possible to find 
experimental situations where the results of these theories are significantly different, and can 
be tested experimentally, is still an open question (for instance, whether or not ISRE in lateral 
collisions may become dominant). 

A second conclusion is that a correct treatment of the correlations between internal vari- 
ables and velocities at each point of space may be important; this is illustrated by the figures 
contained in figure 2. In other words, hydrodynamic equations are not always appropriate. 
What is remarkable is that the effects of local field inhomogeneity and velocities combine non- 
linearly to create this almost complete separation of the two internal states; the effect is not 
a segregation of states but rather an internal conversion that depends on the direction of the 
velocity. It is probably even more remarkable that this effect was discovered experimentally, 
without any theoretical prediction as a guide to the appropriate experimental conditions, and 
appeared basically as a 100% effect from the beginning. Clearly, ISRE is not a small quantum 
correction to the effect of collisions, but may dominate the entire dynamics of the quantum 
gas! 

It would be interesting to extend the present work in a few directions, for instance exploring 
the influence of a temperature gradient on the spin current in hydro dynamical situations, as 
predicted in For bosons, it would also be important to understand better how spin waves 
evolve progressively from the non-condensed regime, which we have studied in this article, to 
the condensed low temperature regime; the theory recently developed in [22] seems to provide 
an appropriate tool for this purpose. 
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APPENDICES: INTERACTION TERM OF THE S-MATRIX ANSATZ 

Appendix I: spin-independent interactions 

In this appendix, we give the detailed calculation of the Wigner transform of the collision 
term of the S matrix Ansatz, equation H18|l. We will closely follow the notations of a similar 
calculation done in the Appendix of Ref. |14j , where the Wigner transform of the Yvon-Snider 
Ansatz, equation @ without internal levels and without exchange, is obtained to first order 
in the gradient expansion. We define the on-shell T matrix T by: 

S = 1- ilnf (58) 

and rewrite the core of the Ansatz in equation ()18() as: 

^pi(l)pi(2)5t - /5i(l)/)i(2) = -i2Trfpiil)pi{2) + c.c. 

+ (27r)2tpi(l)pi(2)tt (59) 
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where c.c. is the complex conjugate of the preceding term. The first two terms of the preceding 
equation are hnear in the T matrix, the last is quadratic. The Wigner transform Ft of the 
linear in T matrix operator between curly brackets in H18|) is: 



Fr(R,r,P,p) 



x(K+,k+|i±^rpi(l)pi(2)i±#i^|K_,k_)+c.c. 



V2 



V2 



where: 



P K , p k 

K± = - ± — ; k± = ^- ± - 

h 2' h 2 



(60) 



(61) 



With two closure relations for relative wavevectors ki and k2, the r.h.s. of H6U() becomes equal 
to: 



d^K / d^K / d^ki / d^k2 e 



iKR„ifcr 



2i{2-K)^h^ 

X (r(k+, ki) + er(-k+, ki)Ps) 5{Ek^ - Ek,] 
x(K+,ki|pi(l)pi(2)|K_,k2) (<5(k2 - k_) + e(5(k2 + \^^)Ps ) + c.c. 



(62) 



where Ps is the exchange operator in spin space. Now using the inverse Wigner transform 
formula, we can introduce into 1)62(1 the Wigner transform /5(r, p) of pi and replace the matrix 
element of the product of /5's by: 



^6 / d?R' J dSr'e-^K-^'e^C^^-kO-r' 



xp 



+ 



P+aki+;ik2 
2 ' 2 



nlTt' — ^ P-^ki-;ik2 
p 2 ' 2 



When inserting the result into ()62j) . a delta function appears: 

I d3i^e^K.(R-R') ^ (27r)3<5(R - R') 



(63) 



(64) 



Tracing over the orbital and the spin space of particle 2, we finally obtain the following 
expression for the Wigner transform It (ri, Pi) for the terms linear in T in the Ansatz H18() : 



-^r(ri,pi 



1 



2i(27r)2 At 
X e """e 



d'q I d-'r I d^K I Sk[ I d'k'^ I d'r' 



73 iJ 



X 5(k'2-k_)r(k+,k;)p(^ri- 
+ ed(k^-k_)r(-k+,k;)/5(^ri 
+ e5(k'2 + k_)r(k+,k;)pfri- 



2 

r + r' 



P2 P ri 



r + r' 
2 ' 
r — r 



P2 



.P'l' 



r — r 



Pi P n 



r + r' 



,P2 



+ 5(k'2+k_)T(-k+,k;)/ n 



Pi /5 ri 



r + r' 



,P2 



C.C. 



(65) 
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with notation (jlll) and with: 



r = ri - r2 ; ^ k± = E ± f ^ 

Pi = Pi - f + ^^V^ ; P2 = Pi - f - ^^T^ 



(66) 



The distribution / is the spin trace of p, see equation (|1U() . 

The same kind of calculation can be done for the term quadratic in the T matrix in p8|). 
It gives the Wigner transform Iy2(ri,pi) as: 



Ir2(ri,pi) 



+ er(k+,k;)r(-k_,k'2)*pfri 



r — r 



r — r 



P2 P ri 



2 

r + r' 



P'l' 



+ r(-k+,k'i)r(™k_,k'2)*/ r 



2 



r — r 



Pi P ri 



r + r' „ 



.P2 



, Pi P ri 



r + r 



,P2 



(67) 



As in Ref. [261 I14j . we now assume that p{r,p) varies slowly in space over microscopic 
distances and expand the product of p and / in and (jHTj) according to 



p(ri,Pi)/(ri,P2) - ^ • [/5(ri,p")/(ri,P2)] 
+ ^ • [p(ri,p'/)VrJ(ri,p'2') - (Vr,/5(ri,p;')) /(ri,p'2') + 



(68) 



The first term in (|68|) corresponds to the local term, the two which follow to first order in 
gradients non-local terms. A very similar equation exists for the product of p's in (|65jl and 

(inzi)- 

(1) Local term We now calculate the zeroth order in gradients. The following three 
integrals occur: 

J d^r' e-^(i'2+'^i) '-' = {2Trf6{k'2 + k'^) (69) 
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The linear in T matrix terms become: 
-^T(ri,pi; 



X 5(k; - k)T(k, k[)p (n, pi - ;ik+;ik;) / (n, pi - hk-hk[) 

+ e5(k; - k)T(-k,k'i)p (ri,pi - hk-nk[) p (ri,pi - nk+hk[) 

+ e6ik[ + k)T(k, k[)p (ri, pi - ;ik+;ik;) p (ri, pi - ;ik-;ik;) 

+ 5{k[ + k)r(-k, k[)f (n, PI - hk+hk[) p (n, PI - hk-nk[) 

+ c.c. 



(70) 



In the preceding equation, we note the appearance of the square of delta functions of the 
energy 6{Ef:)6{k) oc [5{Ek)]'^. They are handled by the following well-known simplification 
from scattering theory: 



[S{Ek 



At 
27rh 



(71) 



where At is a time larger than the duration of a collision (see Ref . ^484 , for example) , which 
simplifies with the one introduced in the S matrix Ansatz (|18)1 . The same manipulations can 
be done with the quadratic in T matrix term, they lead to: 



-^T2(ri,Pi) 



2h 



X r(k,k')T(k,k')*/5(ri,pi -nk + /ik')/(ri,pi -/ik-nk') 

+ eT{-k, k')r(k, k')*p (ri , pi - ^ik - ^k') p (ri, pi - ^k + hk') 
+ eT{k, k')T{-k, k')*p (ri, pi -hk + hk') p (ri, pi - ^k - ^ik') 

+ r(-k, k')r(-k, k')*f (n, pi -nk + hk') p (n, pi -hk- hk') 



(72) 



We will now use different properties of the T matrix. In addition to the differential 
ak{k,k') and total c7t(A;) cross sections defined in (|12() . we introduce, following Ref. pjj, the 
following cross sections: 



r(-k,k) 
r(-k',k)r(k',k)* 



h^k 
h^ 



^rZ:Ak)) 



i'^fwd.i.k) — I'l 

arik,k')-iTr{k,k') 



Unitarity of the S matrix implies the optical theorem: 

r(k,k) 



+ c.c. 



Imr(k,k) 



h^k 



-arik) 



2i ■ (27r)3m 
and the rotational invariance of the interaction Hamiltonian can be used to show that: 

T(-k,-k') =r(k,k') 



(73) 



(74) 



(75) 



Performing the integral over the length of wavevector k' in (|7U() and ()72() and using the 
previous properties of the T matrix, we obtain the Wigner transform to zero order of the S 
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matrix Ansatz Iw = Ij- + It^ ■ 



%(ri,Pi 



m 



cjT(A;)p(ri,pi)/(ri,p2) 



+ o i'^Twd.ik) [pir I, pi), p{r I, p2)], + iT^f^aX^) [p(ri, pi), p(ri, p2)]_ 



ak{k, /c')/5(ri, p'i)/(ri, p^) + - (^af (A:, k') [p(ri, p'^), p(ri, p^) 



+ iTl''-{k,k') [/5(ri,p'i),p(ri,p^) 



(76) 



with the notations This local term is the LL "collision integral" ^Hl- It contains terms 

linear in T matrix, usually written on the l.h.s. of the kinetic equation, as well as terms 
quadratic in T matrix. 

The low-energy limit of the cross sections is obtained from the low-energy expression of 
the T matrix 

(TT(fc) ~ 47ra^ 'yk{k, k') ~ a? 
(^) ~ 47ra2 al'^-{k,k')r^a^ 
rf^Ak) - -^-rra/k rl-\kk')^Q 

when A; — > 0, see Ref. When introduced into (|76j) . the collision integral of equation 

(|22() is obtained, as well as the spin mean field contained in the commutator in equation 
(|2fl|) . Other mean field terms (contained in the anticommutator) appear to first order of the 
gradient expansion. 



(2) First order terms We retain only first order terms that are in addition linear in T 
matrix. 



(2-a) r gradients The first-order terms introduce the gradient of a delta function: 



which implies taking the derivative with respect to k of the function under the integral. We 
obtain: 



(2vr)= 



j A|ReT(k,k) Vp-Vr [p(ri,pi)/(ri,p2)] 

- Vp [Rer(k,k)] • Vr [p(ri,pi)/(ri,p2)] 
+ |r(k,k) Vp • Vr [/5(ri,pi)p(ri,p2)] 

- |Vp [r(k, k)] • Vr [p(ri, pi)p(ri, P2)] 
+ ^r(k, k)* Vp • Vr [p(ri, p2)p(ri. Pi)] 



-Vp [r(k,k)*].Vr [p(ri,p2)p(ri,pi)] 



(77) 
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(2-b) r' gradients The following integrals occur: 

j d?r' e-^(i^2+ki) r'r' = -i(27r)3Vk'/(k^ + k'^) (78) 

so that we obtain: 

^ j A Vp-|ReT(k,k) [p(ri,pi)Vr/(ri,p2)-/(ri,p2)Vrp(ri,pi)] 
+ -k) [/5(ri, pi)Vr/5(ri, P2) - /5(ri, p2)Vrp(ri, pi)] 

+ |T(k, -k)* [Vr/>(ri, P2) /5(ri, pi) - Vr/)(ri, pi) /)(ri, P2)] | (79) 

The expression of the T matrix to lowest order in the scattering length (see equation 0) 

(27r)3r(k,k)~5 

is independent of k. Using the preceding expression in (|77|) and (|79|) . we obtain the following 
first order term: ^ 

2 [Vpp(ri,pi),-Vr(fifn(ri)i + egn{Ti))]^ 

This is the mean field term appearing in the anticommutator in the kinetic equation H2U|) . 

Appendix II: spin-dependent interactions (forward scattering) 

In this appendix, we show how to obtain the full mean field in the case of spin-dependent 
interactions, using the S'-matrix Ansatz. We will limit ourselves to collisions at low-energy, 
using the T matrix only at lowest order in the scattering lengths (equation (|23|) '): 



{a ; P\ fk 17;'^)= ^0^^13,6 + (1 - ^a,l3)Sa,sS, 



3 I V- "a,P/"a,<'"/3,7|-2^-)3 



(80) 



where g^a = Qaa- In this limit the T matrix elements are real. 

As we are interested in mean field terms, we only keep terms linear in the T matrix in the 
r.h.s. of the S matrix Ansatz, equation (|18|) : 



2tt 11 + ePex. 



7 pxiS) Pi(2)-pi(l) />i(2) t 



1 + ^Pe. 

V2 



(81) 



where we introduced the T matrix (the on-shell T matrix), whose definition is given in 
Appendix I, equation Using the properties of the trace on particle 2 and the fact that 

the exchange operator commutes with the T matrix, we can rewrite the preceding equation 
in the form of a commutator: 



1 

ih 



2'Kh 

'At 



rr2{(i + ePex.)tpi(2)}, pi(l) 



(82) 
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It seems therefore natural that the first operator in the commutator should play the role of 
an effective single-particle Hamiltonian. 

We now consider a general matrix element of the preceding commutator. To simplify the 
notation, we only write the first term of the commutator: 

(1 :pi,a|^rr2{(i + ePex.)t /3i(2)}pi(l)|l :p;,a') (83) 

We calculate this matrix element by introducing three closure relations of the form: 

1= / d3p^|p,/5)(p,/3| (84) 

and by using the expression (|8fl|) of the T matrix. We obtain: 

Wf^t J "^'^^ J "^'^^ '^(^(Pl-P2)/2 - '^(p3-P4)/2) 
X Yl \ 5'aj(P3'«IPl|Pl'«')(P4,/3|/5l|P2,/3) 

+ (1 - ^a,(s)9al{P3,P\pi\Pi,a'){P4.,a\pi\p2,f3) 
+ (^9al (P3, Ip'i, a')(p4, a|pi|P2, /5) 

+ e{l - 6a,i3)gal{P3,a\pi\p'i,a'){p4, 13\pi\p2, P) \ (85) 



where P4 = Pi + P2 — Pa- By defining the three following coupling constants: 

yl3 ^ ^9a,f3 



9aa=g^^l ; = + e<?i*lj if « / /? (86) 



we can rewrite equation (|85|) as: 

2TTh 



At 



X 



j d^P2 j d^p3 -5(^(p^_p2)/2 - -E(p3_p4)/2) 

^ (f^|(P4'/5|Pi|P2,/3)(P3,a|/5i|p'i,a'; 



/3 

+ e(P4,a|/5i|p2,/3)(P3,/5|/5i|p'i,a')| (87) 

We therefore see that, for both statistics (bosonic or fermionic), only three coupling constants 
(511 ) 922 and (712) are involved. 

Now, taking the Wigner transform of the operator (|82|) and expanding in gradients (see 
Appendix I) introduces another delta function of energy conservation, which is present at each 
order in gradients. This allows to use the "^(E')^" simplification (equation (|7H) in Appendix 
I) to make the At disappear. The net result is that: 

— (5(^(pi_P2)/2 - -E^(p3-p4)/2) (88) 
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gets replaced by 1. Formally the matrix element of equation becomes: 
J d^p2 j <^^P-i^-^^ X |(P4,/3|/Oi|P2,/3)(P3,a|/5i|p'i,a'; 

+ e(P4, a\pi |P2, /3)(P3, /?l/5i |p'i, a') 
This result can be rewritten: 

(l:pi,a|y-^(l)pi(l)|l:p'i,a') (90) 

where: 

yai3{'^) = ^ap ^ ga^(r,7|pi|r,7) + ec/a/3(r,a|/5i|r,/3) (91) 

7=1,2 

is the mean field potential. Therefore, equation (|82|) can be written formally: 

1[f"^(1), pi(1)]_ (92) 

This equation is valid provided it is used to compute the Wigner transform and do a subse- 
quent gradient expansion. 

In the body of the article, we define the effective single-particle Hamiltonian U which is 
the sum of the external potential y^^* and of the mean field potential V'^^ . The operatorial 
kinetic equation is now: 

^ + {ih)-' + =0 (93) 

at L J - 

The kinetic equation is then derived following the procedure of Appendix I (Wigner transform 
and gradient expansion up to first order); the result is: 

dtpir, P) + - • Vr/5(r, p) + 4[p(r, p), ^(r)]_ 
m in 

-^[Vpp(r,p),.VrC/(r)]+ = (94) 
which is equal to the l.h.s. of ()25|) . 
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